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^ ■ Abstract 

(-1 I We have studied the production of the ^5*0 charmonium state, r]c, at the Large 

Qh' Hadron Cohider (LHC) in the framework of Non-Relativistic Quantum Chromo- 

o ■ dynamics (NRQCD) using heavy-quark symmetry. We find that NRQCD predicts 

O I a large production cross-section for this resonance at the LHC even after taking 

account the small branching ratio of rjc into two photons. We show that it will be 
possible to test NRQCD through its predictions for ijc, with the statistics that will 
be achieved at the early stage of the LHC, running at a center of mass energy of 7 
C^ ' TeV with an integrated luminosity of 100 pb~^. 

Non-Relativistic Quantum Chromodynamics (NRQCD) [I] is an effective theory that has 
been extensively used to study the production and decay of quarkonia. NRQCD is derived 



> 

m 
in 



O I from the QCD Lagrangian by neglecting all states of momenta much larger than the heavy 

^ ■ quark mass, Mq and to account for this exclusion by adding new interaction terms in the 

effective Lagrangian. It is then possible to expand the quarkonium state in-terms of v, 

the relative velocity of the heavy quarks in the bound state. In this expansion, the QQ 

r> I pair in the intermediate state can be in either colour-singlet or colour-octet configurations 

c^ I denoted by QQ[^^~^^Lj ]. However, the colour-octet QQ state evolves non-perturbatively 

into a physical colour-singlet state by emission of one or more soft gluons. The cross 

section for production of a quarkonium state H can be factorised as: 

-iH) = E -v/i:^Ao^C'^'Lj)), (1) 

n={a,S,L,J} ^^-'Q 

where F^s are the short-distance coefficients and On are operators of naive dimension dn, 
describing the long-distance effects. These non-perturbative matrix elements are guaran- 
teed to be energy-independent due to the NRQCD factorization formula, so that they may 
be extracted at a given energy and used to predict quarkonium cross-sections at other 
energies. 
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Before the effective theory approach of NRQCD was developed, the colour-singlet 
model (CSM) [21 [3] was used to analyze the production of quarkonia, where the QQ state 
produced in the short-distance process was assumed to be a colour- singlet. However, it was 
pointed out in Ref. |lj that contributions from the colour-octet operators are significant in 
describing the phenomenology of large-p^ P-state charmonium production at the Tevatron 
[5]. In Refs. [6], [7] the complete set of short-distance coefficients in NRQCD needed to 
study J /ip and x production was calculated and compared with the data from Tevatron lj. 
These NRQCD calculations gave a good description of the shapes of the p^ distributions of 
the charmonium resonances at the Tevatron but the normalization of these distributions 
was not predicted in NRQCD i.e. the non-perturbative matrix elements which determined 
the normalization had to be obtained by a fit to the data. Independent tests of the effective 
theory approach were, therefore, necessary to determine the validity of the approach and, 
indeed, various proposals were made [TO] to test NRQCD. But several of these proposals 
are not for large-p^, quarkonium production and the validity of NRQCD factorization at 
low-py is suspect. 

One interesting test of NRQCD comes from the study of the polarization of J/^'s at 
large-pj, [11] which primarily comes from a fragmentation-like processes where a single 
gluon splits into a QQ pair which inherits the transverse polarization of the gluon. The 
heavy-quark symmetry of NRQCD then comes into play in protecting this transverse 
polarization in the non-perturbative evolution of the QQ pair into a J /ip. The large-p^ 
J /ip is, therefore, strongly transversely polarized. This is not true at even moderately low 
Py where the J /ip is essentially unpolarized. The p^ dependence of the polarization is, 
therefore, a very good test of the theory [12]. 

The CDF experiment has measured the p^-dependence of the polarization and they 
find no evidence for any transverse polarization at large p^ [13], which seems to indicate a 
dramatic failure of the theory. Inspite of the successful prediction of the production cross- 
sections of the various charmonium resonances it may well be that the effective theory 
is missing out on some aspect of the physics of quarkonium formation. Alternatively, it 
could be that the charm quark is too light to be treated in NRQCD. On the other hand, 
polarization measurements are usually fraught with problems and it may well be that 
the problem is elsewhere. Finally, because the colour-singlet channel predicts unpolarized 
J/V^'s, there have been attempts to increase up the colour-singlet contribution to the 
production processes by invoking Reggeized gluons [H] or enhanced effects of higher- 
order QCD corrections in the singlet channel [151 [IS]- For reviews of the current status 
of these calculations and their experimental consequences, see Refs. [TTlfTS] . 

In this situation, it is worthwhile looking for other tests of NRQCD which success- 
fully navigate between low-p^ and polarization. The heavy-quark symmetry of NRQCD 
provides a set of relations between non-perturbative parameters of different resonances 
so a measurement of a given state yields information on the non-perturbative parameter 
of another state related to the former by heavy-quark symmetry. This fact has been 
exploited to study h^ production at the Tevatron [19j and, more recently, at the LHC 
|2U] . Similarly r^^ production at the Tevatron has also been studied [21j. In this paper, 
we study the production of r]c at the LHC. 



"^See also Ref. [8 . For a detailed review of quarkonium production see Ref. [S] 



The Fock space expansion of the physical ?7c, which is a ^Sq ( J^'" = +) state, is: 
\Vc) = 0(1) gg[^4'']) + Oiv') QQl'Pl') g) + Oiv"") QQl'S?] g) + ■ ■ ■ . (2) 



In the above expansion the colour-singlet ^5*0 state contributes at (9(1). As the P-state 
production is itself down by factor of (9(f^) both the colour-octet ^Pi and ^Si channels 
effectively contribute at the same order. The colour-octet state ^P| (^>S'| ) becomes a 
physical rjc by emitting a gluon in an El (Ml) transition. Keeping terms up-to (9(aff ^) 
the r]c production cross section can be parameterized as: 

c^iVc) = ^^ {0\OnSo]\0) 

+^^ (0| O^Pi] |0) + ^^ (0| O^fS^] |0) , (3) 

where the coefficients, F's, are the cross sections for the production of cc pair in the 
respective angular momentum and colour states. The differential cross section for cc pair 
production with specific angular momentum and colour states at the LHC is given by: 

^ {pp ^ c-c p^+^L[^'«]] X) = 

Y, dy dxi xi Ga/p{xi) X2 Gb/p{x2) h — - -jji(^b -^ cc[^^+^l5'^'] d), (4) 

■J J Zx\ — Xf c^ dt 

where the summation is over the partons (a and b), the final state QQ is in the ^Sq , 
^Pi , ^Si states and Ga/p, Gb/p are the distributions of partons a and b in the protons 
and Xi, X2 are the respective momentum they carry. X2 is related with Xi as: 

_ xi XT e-y - 2t 
^^ ~ 2xi-XTey ' ^ ^ 



where xt = Jxj^ + 4r = 2Mt/ \fs with xt = 2p^j ^ and r = M"^ / s. Here a/s is the 
center-of-mass energy, M is the mass of the resonance and y is the rapidity at which the 
resonance is produced. The subprocesses contributing to Eq.(|l]) are: 

gg ^ QW'^'Lj]g, 

gq{q) ^ QW'^'Lj] q{q), (6) 

qq ^ QQ[^'^'Lj]g. 

The matrix elements for the subprocesses corresponding to Fi[^S'o] and FspS*!] are listed 
in Refs. [71|22]. The remaining coefficient F8[^Pi] has been calculated and used in [21] to 
analyze rjc production at the Tevatron. 

We use heavy quark spin-symmetry to obtain the values of {O^Ys from the experi- 
mentally predicted values (C;(/'^)'s. Using this symmetry the {0^)^s are related as: 

{0\Or['So]\0) = l{0\Oi^^['S^]\0) {l + Oiv')l 

(0|OrrPi]|0) = (0|08^/^pPo]|0) (l + 0(^;2)), 

(0|OrP^i]|0) = {0\Oi^^['So]\0) {l + 0{v')). (7) 



We use the predicted value of the singlet matrix elements listed in Ref. |6] ( Ci = 
(0|C/'^'^pS'i] |0) = 1.2 GeV^) and the octet matrix elements extracted from the CDF 
J/i^ data I7j (i. e. C, + C,^ mt/l^m + fM^^l^m. = (2.2 ± 0.5) x IQ-^ GeV^). It is 



to be noted that only a linear combination C2 + C3 can be extracted from the CDF data 
as the shapes of C2 and C3 contributions to the J/ip p^-distribution are almost identical. 
Hence, for predicting the rjc rate we assume that either C2 or C3 saturates the sum so 
that our predictions indicate the band within which we expect the experimental value of 
the 7]c production cross-section to lie. More explicitly, we consider the two extreme cases 
where in the first case, the maximum possible contribution is from the '^5'| channel and 



none from the ^P| channel whereas in the second case the ^P| contributes its maximum 
while the ^S[ channel does not contribute. 
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Figure 1: da/dpj, (in nb/GeV) for rjc production (after folding in with Br(?7c — ^ 77) = 
3.0 X IQ-^) in pp collisions at v^i = 7 TeV and 14 TeV with -2 < y < 2. 

In Fig. [1] we have displayed the differential cross section Bda/dp^ as a function of Pj, 
at two different center-of-mass energies, viz. -y/i = 7 and 14 TeV respectively, where B 
is the r]c — )■ 77 branching ratio (i? = 3 x 10~^). We have used CTEQ 5L LO parton 
densities [23] evolved to a scale Q = Mt- In both (a) and (b) of Fig. 1, we have curves 
marked I and II where I is the sum of the colour- singlet and the ^S*} contributions and II 
is, likewise, the sum of the colour- singlet and the ^P^ contribution. In addition, we also 
display the colour- singlet curve in both the figures, to bring out the fact that the octet 
contributions overwhelmingly dominate the cross-section. 

To account for the experimental threshold in p^ for the photons which is about 10 
GeV, we use a lower-p^, cut of 20 GeV on the rjc to calculate the integrated cross-sections. 
For the LHC running at ^/s = 7 TeV, with an integrated luminosity (£) of 100 pb~^, we 
find that the number of events in the 77 channel from the singlet ^Sq is about 40 while 
the number of events from ^5*1 is about 10^ when ^P^ contribution is neglected, while 
the number of events from ^P| is about 37660 when ^5*} contribution is absent. For the 



case of y/s = 14 TeV and C = 100 pb ^, the respective numbers would be 100, 2.8 x 10^ 



and 10^. Thus the minimum rjc — ?■ 77 events at the LHC running at y/s 
will be 37700 and 10^ respectively. 



7 and 14 TeV 



We can see from Fig. [T]that the shapes of the ^S[ and the ^P^ contributions to the 
Pj,-distribution is different and may allow the non-perturbative matrix elements C2 and 
C3 to be individually determined. We have noted earlier that the J/ip cross-section does 
not provide such a separation. But it is also pertinent to note that the integrated cross- 
sections are also very sensitive to the the values of C2 and C3 and so the measurement of 
the integrated cross-section alone may provide this discriminatory ability. 
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Figure 2: Variation of total cross section with respect to chosen minimum p^-cut for rjc 
production (after folding in with Br(r7c -^ 77) = 3.0 x 10"'^) in pp collisions at y/s = 
7 TeV and 14 TeV with -2 < y < 2. 



In Fig. [2] we have displayed the effect of increasing the Pr^-cnt on the magnitude of the 
integrated cross-section. As explained earlier, the cut on p^ of the 7]^ will be determined 
by the minimum p^ threshold that the experiments use to trigger on the photons, for 
which the usual choice is 10 GeV. In case the experiments use a larger cut on the p^ of 
the photons in order to improve the quality of their signal, the p^, cut on the rjc will be 
correspondingly higher. We see, from Fig. 2, that even when the cut on the minimum p^, 
is as large as 50 GeV, the cross-section is substantial. 

We also have analyzed the effect on the cross-section of the variation of the QCD 
scale, the parton densities and the non-perturbative matrix elements. Table [1] shows the 
variation in the minimum and maximum number of 7]c events with scale Q, expected at 
the LHC running at a/s = 7 and 14 TeV. We find that the cross section decreases by 
about 25-30% changing the scale from Q = Mt to Q = 2Mt and it increases by 40-50% 
for the scale choice Q = Mt/2 instead oi Q = Mt- We have checked that the cross section 
decreases by about 10-20% if we use MRST LO densities [21] instead of CTEQ5L LO [23j 
densities. Since the heavy-quark symmetry is an approximate symmetry we can expect 
about 30% variation in the values of the non-perturbative matrix elements we have used. 





~ Number of rjc events evaluated using CTEQ 5L LO densities 


v^ 


Q = Mt/2 Q = Mt Q = 2 Mt 


7TeV 


5.6 X 10^ - 1.5 X 10^ 3.8 x 10^ -10^ 2.6 x 10^ - 7.3 x 10^ 


14TeV 


1.5 X 10^ - 3.9 X 10^ 10^ - 2.8 x 10^ 7.6 x 10^ - 2.1 x 10^ 




~ Number of events evaluated using MRST LO densities 


v^ 


Q = Mt/2 Q = Mt Q = 2 Mt 


7TeV 


4.9 X 10^ - 1.3 X 10^ 3.4 x 10^ - 9.5 x 10^ 2.4 x 10^ - 6.8 x 10^ 


14TeV 


1.2 X 10^ - 3.3 X 10^ 9 X 10^ - 2.4 x 10^ 6.8 x 10^ - 1.8 x 10^ 



Table 1: Minimum and maximum number of rjc events expected at the LHC for an 
integrated luminosity of 100 pb~^. 



All through we have considered only the direct production of rjc at the LHC. However, 
an additional contribution to r]c signal comes from the decays of J/ip. We have estimated 
that this additional contribution to the signal coming from J/ip decays can change our 
predictions by about 1% as Br(J/-?/' — )■ rjc'y) ~ O{10~^) and the J/ip production cross 
section is expected to be of same order of r]c production cross section. 

We would like to remark that such an analysis may also be carried out for the bot- 
tomonium resonance r^f,. The corresponding non-perturbative parameters in that case, 
however, are very poorly determined and suffer from large errors. 

In conclusion, the heavy-quark symmetry of NRQCD allows us to make predictions 
for rjc production at the LHC. Measurements of the integrated cross-section and the p^, 
distribution of rjc at the LHC will provide a very good test of NRQCD. We show that 
NRQCD predicts a large cross-section for rjc at the LHC even at i/i = 7 TeV and so this 
prediction is easily testable. 



Acknowledgments 

We thank Gobinda Majumder and Gagan Mohanty for useful discussions. K.S. completed 
part of the work presented here on a visit to the Department of Physics and Astronomy, 
University of Southampton, U.K. He gratefully acknowledges support from the University 
of Southampton for this visit. 



References 

[1] G. T. Bodwin, E. Braaten and G. P. Lepage, Phys. Rev. D 51 (1995) 1125 [Erratum- 
ibid. D 55 (1997) 5853] [arXiv:hep-ph/940733 9|. 

[2] E. L. Berger and D. L. Jones, Phys. Rev. D 23 (1981) 1521. 

[3] R. Baier and R. Ruckl, Z. Phys. C 19 (1983) 251. 

[4] E. Braaten, M. A. Doncheski, S. Fleming and M. L. Mangano, Phys. Lett. B 333 



(1994) 548 arXiv:hep-ph/9405407]; D. P. Roy and K. Sridhar, Phys. Lett. B 339 



(1994) 141 |arXiv:hep-ph/9406386|; M. Cacciari and M. Greco, Phys. Rev. Lett. 73 
(1994) 1586 |ar Xiv:hep-ph/9405241i . 

[5] F. Abe et al. [CDF Collaboration], Phys. Rev. Lett. 69 (1992) 3704; Phys. Rev. 
Lett. 79 (1997) 572; Phys. Rev. Lett. 79 (1997) 5 78; D. E. A costa et al. [CDF 
Collaboration], Phys. Rev. D 71 (2005) 032001 [a rXiv:hep-ex/0412071j . 

[6] P. L. Cho and A. K. Leibovich, Phys. Rev. D 53 (1996) 150 |arXiv:hep-ph/95 05329]. 

[7] P. L. Cho and A. K. Leibovich, Phys. Rev. D 53 (1996) 6203 | arXiv:hep-ph/9511315| . 

[8] M. Cacciari, M. Greco, M. L. Mangano and A. Petrelli, Phys. Lett. B 356 (1995) 
553 [arXiv:hep-ph/9505379|. 



[9] N. Brambilla et al. [Quarkonium Working Group], arXiv:hep-ph/0412158 



[10] M. Cacciari and M. Kramer, Phys. Rev. Lett. 76 (1996) 4128 [a rXiv:hep-p h/9601276|; 
J. Amundson, S. Fleming and I. Maksymyk, Phys. Rev. D 56 (1997) 5844 
|arXiv:hep-ph/9601298|; S. Gupta and K. Sridhar, Phys. Rev. D 54 (1996) 5545 
[arXiv:hep-p h/9601349j ; Phys. Rev. D 55 (1997) 2650 |arXiv:hep-p h/960843 3j; 
M. Beneke a nd L Z. ploths tein, Phys. Rev. D 54 (1996) 2005 [Erratum-ibid. D 54 
(1996) 7082] [arXiv:hep-ph /9603400 1 ; W. K. Tang and M. Vanttinen, Phys. Rev. D 
54 (1996) 4349 [arXiv:hep-ph/9603266|; E. Braaten and Y. Q. Chen, Phys. Rev. Lett. 
76 (1996) 730 |arXiv:hep-ph/9508373|; K. M. Cheung, W. Y. Keung and T. C. Yuan, 
Phys. Rev. Lett. 76 (199 6) 877 |arXiv:h ep-ph/9509308 |; P. L. Cho, Phys. Lett. B 368 
(1996) 171 |arXi v:hep-p h/9509 355|; K. M. Cheung, W. Y. Keung and T. C. Yuan, 
Phys. Rev. D 54 (1996) 929 [a rXiv:hep-ph/9602423| ; P. Ko, J. Lee and H. S. Song, 
Phys. Rev. D 53 (1996) 1409 [arXiv:hep-ph/9510202|; G. T. B odwin, E. Braaten, 
T. C. Yuan and G. P. Lepage, Phys. Rev. D 46 (1992) 3703 ^arXiv:hep-ph/9208254j . 



[11] P. L. Cho and M. B. Wise, Phys. Lett. B 346 (1995) 129 arXiv:hep-ph/9411303 



[12] M. Beneke and M. Kramer, Phys. Rev. D 55 (1997) 5269 arXiv:hep-ph/9611218]. 



[13] A. A. Affolder et al. [CDF Collaboration], Phys. Rev. Lett. 85 (2000) 2886 
[arXiv:hep-e x/0004027] ; A. Abulencia et al. [CDF Collaboration], Phys. Rev. Lett. 
99 (2007) 132001 [a rXiv: 0704. 0638 1 [hep-ex]]. 

[14] V. A. Khoze, A. D. Martin, M. G. Ryskin and W. J. Stirling, Eur. Phys. J. C 39 
(2005) 163 [arXiv:hep-ph/0410020 |. 



[15 
[16 

[ir 

[18 
[19 
[20 
[21 

[22 
[23 



B. Gong, X. Q. Li and J. X. Wang, larXiv: 0805 .47511 [hep-ph]. 

P. Artoisenet, J. Campbell, J. P. Lansberg, F. Maltoni and F. Traniontano, Phys. 
Rev. Lett. 101 (2008) 152001 [arXiv:0806.3282 [hep-ph]]. 

J. P. Lansberg et ai, AIP Conf. Proc. 1038 (2008) 15 [ar Xiv:0807.3666l [hep-ph]]. 

J. P. Lansberg. larXiv:081L4 005 [hep-ph]. 



K. Sridhar, Phys. Rev. Lett. 77 (1996) 4880 |arXiv:hep-ph/9609285] . 

K. Sridhar, Phys. Lett. B 674 (2009) 36 [arXiv:0812.0474 f [hep-ph]]. 

P. Mathews, P. Pou lose and K. Sridhar, Phys. Lett. B 438 (1998) 336 
|arXi v:hep-p h/9803424| . 

R. Gastmans, W. Troost and T. T. Wu, Nucl. Phys. B 291 (1987) 731. 

H. L. Lai et ai [GTEQ Collaboration], Eur. Phys. J. C 12, (2000) 375 
|arXiv:hep-ph/9903282]. 



[24] A. D. Martin, W. J. Stirling and R. G. Roberts, Phys. Lett. B 306 (1993) 145 
[Erratum-ibid. B 309 (1993) 492]. 



